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This work investigates how a conical singularity can affect the specific heat of systems. A free 
nonrelativistic particle confined to the lateral surface of a cone — conical box — is taken as a toy 
model. Its specific heat is determined as a function of the deficit angle and the temperature. For 
a vanishing deficit angle, the specific heat is that of a particle in a flat disk where a characteristic 
temperature separates quantum and classical behaviors, as usual. By increasing the deficit angle the 
characteristic temperature increases also, and eventually another characteristic temperature (which 
does not depend on the deficit angle) arises. When the cone gets sufficiently sharp, at low and 
intermediate temperatures the azimuthal degree of freedom is suppressed. At low temperatures the 
. . . specific heat varies discontinuously with the deficit angle. Connections between certain theorems 

>0 ' regarding common zeros of the Bessel functions and this discontinuity are reported. 
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I \ As techniques in condensed matter physics to fabricating very thin samples with various topologies and geometries 

J> ■ develop, the study of quantum mechanics on curved surfaces becomes increasingly relevant. Among these backgrounds 

0^ ' a cone is perhaps the one with the simplest nontrivial geometry — it is fiat everywhere except at its tip where there 

'""' , is a conical singularity [ll Qj S lij ^J • Due to local flatness the conical singularity does not affect local dynamics; but 

^P ' it does affect global features such as the solutions of differential equations. One is dealing, therefore, with a geometric 

analog of the Aharonov-Bohm setup [g ■ 

There has appeared a considerable number of works on quantum mechanics of a particle in conical spaces, since 

f—^ \ the early investigations in Refs. |3) S I3- Nevertheless, it seems that a study on how a conical singularity affects 

the statistical mechanics of a particle has not yet been pursued. It should be mentioned that since physics around 

a topological defect in the bulk of a solid can be seen as physics in an effective conical geometry [l^, the study 
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I. INTRODUCTION 



o 

O^' . . . . . 

I of statistical mechanics in conical spaces may help to learn about thermal properties when topological defects are 

wp. present. 

• • , This paper addresses quantum statistical mechanics of a free nonrelativistic particle on the cone. The focus is on 

. ^ ' ^^^ specific heat C. In Sec. II the classical partition function is calculated, showing that classical C is not affected by 

Jv>( \ the conical singularity. In Sec. Ill, the particle energy spectrum is obtained and an analysis on how it depends on the 

j_j • deficit angle is presented. In particular, connections with properties of the zeros of the Bessel functions are discussed. 

Ci ' In the following section, the energy spectrum is used to obtaining the quantum specific heat. Quantum C is then 

studied at low, intermediate and high temperatures, showing that the conical singularity can affect it substantially. 

In Sec. V the paper is closed with a summary and final remarks. 

II. CLASSICAL STATISTICAL MECHANICS 

One begins by defining the background (a more detailed account on the conical geometry can be found for example 
in Refs. [a, Ell)- A wedge of angle V is removed from a disk of radius a, and the two opposite sides of the missing 
wedge are identified, resulting a cone of deficit angle T). It is convenient to express T) in terms of the cone parameter 
<a< 1, 

X> == 27r(l - a). (1) 

Then a flat disk (i.e., absence of conical singularity) corresponds to a = 1, and the cone becomes increasingly sharper 
as a decreases from unity. It is worth mentioning that a cone with I? = tt corresponds to the configuration space of 
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a system of two identical particles (here separate by a maximum distance a) with respect to its centre of mass [l^ • 

A free norelativistic particle of mass M in thermodynamical equilibrium with a reservoir at temperature T is 
constrained to move on the cone described above. The corresponding Hamiltonian is given by 

2 2 

'W — ^P -X- ^^ (9\ 

where Pp and p^p are the usual momenta associated with polar coordinates p and (p. Although Eq. ^ carries no 
trace of the conical singularity [Eq. Q is simply the Hamiltonian of a free particle on the Euclidean plane] , the latter 
appears in the identification 

(p,(^ + 27ra) - (p,(/?). (3) 

It can be shown that the conical singularity at p = corresponds to an integral curvature proportional to 2? [j| . 

Observing Eo. (1^ . the standard procedure to compute the (semi) classical partition function by integrating over 
the phase space [l3|. 

^(2) = T^ '^P d^ dpp dpp e-^l'^^, (4) 

" JO JO J-oo J-oo 

yields 

Z(2) = ^It^MkT, (5) 

where 

A := a-KO^ (6) 

is the area of the cone. Once a partition function Z is known, the internal energy is given by 

U = nT'^^. (7) 

Z dT ^ ' 

As the dependence of Z(2) on a happens only through a temperature independent factor, the internal energy and, 
consequently, the specific heat 

are clearly not affected by the conical singularity. Indeed, these formulae give 

[/ == kT C = K (9) 

as one might have guessed by evoking the equipartition theorem at an early stage [see Eq. ^]. [In fact, as a 
straightforward covariant calculation can show, the results in Eqs. (O and (Q hold for any two-dimensional curved 
box regardless its topology and geometry] 

III. ENERGY SPECTRUM 

In order to use the tools of quantum statistical mechanics, one needs to know the energy spectrum of the system 
— a free particle in a conical box. Considering Eq. (j2l the stationary states ^ must satisfy 

^{p,Lp + 2TTa) = *(p, (^). 

Eq. (O leads to a Hamiltonian operator which is simply — ft^/2M times the flat Laplacian written in polar coordinates. 
One then faces the task of solving the eigenvalue problem for the flat Laplacian operator on the cone ^ J] , with the 
additional boundary condition 

*(p = a,(p) = 0, 



since the particle is confined to a circular box of radius a. (This is essentially the eigenvalue problem which appears 
in the study of vibrations of a circular membrane |14|.) 

Up to a normalization constant, the regular solutions (i.e., finite at p = 0) are 

"^ = J\n\/a{kn.sP)exp{in(p/a} kn,s ■■= -X\n\/a.s (10) 

a 

where n is an i nteger s is a positive integer, and Xv,s denotes the sth positive zero of the Bessel function of the first 
kind Jy (x) |l4L llSl Il6l Il7| . Approximated values for the zeros can be obtained by truncating the expansion |l4| 

i^ 1\ -iv^ -I fv 1^ "^ 



at some point which depends on the accuracy required. Eq. (|11() is specially adequate in dealing with high zeros 
(large s). 

Then, for a pair of quantum numbers (n, s), ^ is a state with angular momentum nh/ a (n > and n < correspond 
to the particle winding counterclockwise and clockwise, respectively) and energy 

E^^^t^, (12) 

2M ' ^ ' 

both quantities, therefore, discrete. As En,s = E-n.s, En^Q,s is at least doubly degenerate (degeneracy will be 
addressed in more detail below). 

In order to have a picture of the energy spectrum, and of the corresponding stationary states, one recalls that 
Ju{x) is an oscillatory function of x and that its zeros Xu.s are labelled such that Xv.i < Xv,2 < Xv.z < ■ ■ ■■ It is also 
important to know that Xi^,s is a continuous increasing function of v 15J. Then En^s ^ ^o,1j i-6-j 

ft2 /2.4048...\^ , , 

^°^^ = 2M [-^^) (1^) 

is the zero point energy (xo,i has been replaced by its numerical value jig), and En^s increases with \n\/a and s. 
Clearly, energy levels with n = are not affected by the conical singularity (even classically the conical singularity 
only affects the motion if the particle has nonvanishing angular momentum 11]). 

To figure out the way that a [T>, cf. Eq.lQJ] shapes the energy spectrum, it helps to know that for a — 1 (particle 
confined to a flat disk [13) the energy spectrum resembles that of a particle in a square box of the same area, although 
this analogy should not be pushed too far (e.g., the degeneracies of the corresponding energy levels do not match 
completely). When u > 0, Xu,s > v [13, leading to 

Since E\^\ is the energy of the first excited state with n ^ 0, as the cone gets sharper (i.e., as a decreases from unity) 
all the energy levels with nonvanishing angular momentum are swept up in the spectrum, leaving those with vanishing 
angular momentum unchanged. It should be mentioned that for large v^ Xv,\ — v ^ 1.8557;^"'^/^ [ig, resulting that 
">" in Eq. (|14|) can be replaced by "~" when a is small, with increasing accuracy as a ^ 0. 

In order to discuss degeneracy of the energy levels, it is appropriate to begin with a = 1. As two Bessel functions 
of the first kind with integer orders cannot have common zeros (this follows from a theorem known as Bourget's 
hypothesis [ia|)j it results that energy levels with n = are not degenerate, whereas those with n 7^ are doubly 
degenerate. This situation changes when the conical singularity is present {a ^ 1}, leading to a new source of 
degeneracy, since the Bessel functions in Eq. I|10ll may have common zeros (see Ref. [l^ and references therein). 

Noting that xo,i < Xi,i < X2,i < Xo,2 < . . . 16], when a = 1 Eq. ifT^ yields 

£^0,1 < £^1,1 < £^2,1 < £^0,2 < ■ ■ • ■ (15) 

As a decreases, E\^\ and E2,\ increase, whereas £^0,1 and £^0,2 remain the same. Eventually i?2,i catches £^0,2 up, and 
the energy level £^0,2 — £■2,1 becomes triply degenerate. By further lowering a, the configuration in Eq. (|15|l evolves 
to 

-£0,1 < -£ia < -£0,2 < -£24 < . . . . (16) 



Carrying on with the process, as a reaches the value (numerically determined) 

ao = 0.4338 . . . , 

Eq. (fT?^ changes to 

£'0,1 < £^1,1 = -E'0,2 < £^2,1 < • • ■ , (17) 

where the energy level i?o.2 — Ei^i is again triply degenerate. When a < ao, Eq. H17() is replaced by 

Eq,i <Eo,2<...< Eij <.... (18) 

According to Eq. 114|l and the discussion that follows it, as a ^ more and more levels Eq^s appear on the left of 
Ell in Eq. (|18|l . resulting in that the lower part of the energy spectrum corresponds to states with vanishing angular 
momentum. 

These considerations are relevant to finding the behavior of the specific heat, as will be seen in the next section. 

IV. SPECIFIC HEAT 

In quantum statistical mechanics .13] the internal energy is given by Eq. with the partition function now 
obtained by summing over the stationary states discussed in the previous section, 

00 00 

n— — 00 s— 1 

As usually happens, a closed expression (in terms of known analytic functions) is not available for the summations in 
Eq. p9(l . and a combination of analytical and numerical procedures is required to deal with them. It is convenient 
to define two characteristic temperatures, namely: 

e„ := Ehl^^hl e := ^^bl^3hl (20) 

It follows from the discussion in Sec. Ill that 8q. increases as a decreases from unity, Qa=aa — ©i ^-^d as a ^ 0, 

0"^TTIF^ 2- (21) 

As seen in the previous section, the degree of degeneracy of the first excited energy level is 2, 3 or 1, for a > ao, 
a — ao and a < ao, respectively. Taking this fact to Eq. H19|l . Eqs. lO and lO yield for the low temperature behavior 
(T -^ 0) of the specific heat C{a), 

C{a>ao)^2J^Ve-^-/^ C{a = ao) o, sJ^Y e-''/^ C(a < a„) c. «: (^l V e'^/^ (22) 

where the integers multiplying n are the degree of degeneracy just mentioned. Eq. H22I) shows that C varies discon- 
tinuously with a when T is much smaller than the characteristic temperatures in Eq. (|20l) . As T increases, higher 
energy levels have to be considered in Eq. H19I) . and Eq. H22I) no longer holds. 

The way that C departs from Eq. (|22|l is shown in Fig. ^ In particular, for small values of a the ways that C 
departs from C{a < ao) in Eq. H22|) are approximately the same. This can be understood by recalling that when a is 
small the lower part of the spectrum is essentially formed by energy levels -Eo,s which do not depend on a. Eq. H21|) 
shows that for small a the azimuthal characteristic temperature 0q is high, and consequently the various plots of C 
split from each other only at higher temperatures, when the populations of energy levels with nonvanishing angular 
momentum become appreciable. 

At the regime T <^ 0^ with a < ao [a > ao necessarily leads to the corresponding low temperature behavior in 
Eq. if^ ]. only states with n = effectively counts in Eq. ((T^ . resulting 

00 
Z^J2exp{-Eo,s/KT}. (23) 

s=l 



By setting i^ = in Eq. (HTJ, Eq. lO yields 

which up to the term 1/4 and the corrections is the energy of a free particle in a one-dimensional box of length a 
(and this identification improves as s increases). As T — > 0, Eq. (|23() leads to C{a < ao) in Eq. (|22(l . When a is 
sufficiently small such that 

e < T < e^ (25) 

applies, the summation in Eq. H23I) can be replaced, in first approximation, by an integration over s from 1 to oo, 
with Eq^s given by Eq. 124|) [one can safely truncate Eq. (|24|l at the first term to obtain the main contribution]. 
Using for example Refs. [l^ or [23|, evaluation of the integration yields Z ~ -^(i), where 

Zn) = -^2-kMhT (26) 

^ ' h 

is the classical partition function of a particle in a one-dimensional box of length a. Eqs. and then give 
U ~ kT/2 and 

C-f, (27) 

showing that the conical singularity indeed suppresses one degree of freedom [cf. Eq. ©]. 

With a smaller than about 0.1, Fig. ^ shows that C gets a maximum value just above k/2, at T close to Q (see 
Fig. ^for more accurate values). If 8 ^ 8^ also holds, the regime T <C 8q, resembles the thermal behavior of the 
specific heat of a particle in a one-dimensional box of length a 21]. 

When the temperature is much larger than both the characteristic temperatures, one can proceed as above and 
replace the summations in Eq. (|19|l by integrations. Noting Eqs. (|10|) to H12|) . and performing the integration over s 
|l6l |20j leads to an integration over n of a complementary error function, namely 

I erfc(a;) dx = —=, 
Jo V"' 

which has been evaluated by regularizing related formulae in Refs. [I6l ri7| . The leading contribution in Z is 
proportional to the classical partition function in Eq. ((SJ , with the proportionality constant depending on where Eq. 
(fTT)l is truncated: keeping only the first term, and then also the second term, results 

Z ~ 0.81 Z(2) Z ~ 0.94 Z(2) (28) 

respectively (in fact, to obtain the second expression one has to evaluate an additional integration involving the 
complementary error function [23). Using either results in Eq. (|28|l . Eqs. iQ and © yield as leading contribution 

C~K (29) 

which is the classical value, as expected. [The discrepancies in Eq. (|28() and its failure in reproducing Z^2) exactly 
can be understood by noticing that for large n in Eq. (|19|) . terms beyond the first term in Eq. Hill) become relevant.] 
A maximum also appears, now above k, at a temperature comparable to 8^, as shown in Fig. ^ In particular the 
plot for a = 1 (particle in a flat disk) resembles that of a particle in a square box of the same area |23 . 

V. CONCLUSION 

Motivated by the fact that statistical mechanics of systems in atomic size boxes is very sensitive to the topology 
and geometry of the boxes, this work treated the specific heat C of a free nonrelativistic particle in a conical box. The 
energy spectrum was determined and its features were obtained from properties of the zeros of the Bessel functions. 
In particular, degeneracy of the energy levels was related to long established theorems regarding common zeros of the 
Bessel functions. 

This toy model exposed with clarity the effects of a conical singularity on C at various regimes of temperature. At 
low temperatures C varies discontinuously with the deficit angle [cf. Eq. (|22|l ]. At low and intermediate temperatures, 




T/0 



FIG. 1: Numerical plots of C/k versus T /Q. The plots from the left correspond to a=l, 0.3, 0.2, 0.15,0.1 and a —* Q, respectively. 
The plots for a—1 to a=0.15 display maxima corresponding to specific heat of about 1.119 k, 1.125 k, 1.068 «; and 1.042 n at 
temperatures of about 0.52 0, 0.72 O, 1.59 Q and 3.49 O, respectively. At T ~ 0.94 0, the plot for a -^ displays a maximum 
corresponding to C ~ 0.535 k. 

if the cone is sufRciently sharp, the azimuthal degree of freedom is suppressed [cf. Eq. (|27|l ]. At high temperatures the 
radial and the azimuthal degrees of freedom are fully excited, leading to the value of C prescribed by the equipartition 
theorem [cf. Eq. H29(l ]. The numerical plots in Fig. ^ illustrate these facts. 

The effect of the sharpness of the cone on C can be better appreciated by contrasting the two limits a ~ 1 and 
a ^ 0. When a — I, the energy spectrum is that of a free particle in a flat disk, where energy levels with vanishing 
and nonvanishing angular momentum are mixed up along the spectrum. As temperature T rises, the radial and the 
azimuthal degrees of freedom get equally excited, leading to the high temperature value C ~ k, which is consistent 
with the equipartition theorem. This picture changes radically as one sets a — > 0. When the limit a ^ is literally 
taken, the energy of the first excited state with nonvanishing angular momentum gets infinitely high in the spectrum. 
Then, as T increases, only the radial degree of freedom gets progressively excited, resulting in C ~ k/2 at higher 
temperatures — the a ^ limit suppresses the azimuthal degree of freedom. 

This work also aimed to tune mathematical tools in quantum statistical mechanics on the cone such that other 
thermodynamical quantities, as well as more elaborate models, can be addressed. Some remarks are in order. 

The regular solutions in Eq. IIIUII are not the only ones consistent with square integrability of the wave function 
and conservation of probability 1221 1231 . A mild divergence at the origin (where the conical singularity sits) does not 
spoil these physical requirements [lll | . It seems that the choice for regular solutions relies in that they arise naturally 
when the conical singularity is "smoothed" by regularization |22l |. 

The discontinuity of the specific heat in Eq. H22|) is a feature of any thermal system whose degree of degeneracy 
g of the first excited energy level ei varies with some parameter. If eo is the zero point energy, as T — + 0, a quick 
calculation shows that 



C^gn[- 



-e/T 



where 9 := (ei — eo)/ n- This kind of discontinuity at low temperatures also appears, for example, in the specific heat 
of the anisotropic rigid rotator j2Jj . 

Recalling that even the main contributions in Eq. (|28|l were affected by where one truncates Eq. (|ll|l . the 
determination of corrections in Eqs. H27|l and (|29|) turns out to be a rather tricky task. Although the leading 
contributions in Eqs. (|27|l and (|29l) do not depend on where Eq. (|11|1 is truncated (either on the lower limit of the 
related integrations), this may not be the case when dealing with their corrections. Such a problem is somewhat 
similar to that where two distinct approximations for the energy spectrum give the same leading contribution in the 
expression for the high temperature behavior of C; but different corrections (e.g., the anharmonic oscillator in Ref. 
[25|). Therefore, in order to check correctness at every step in using Eq. Hll|l to finding subleading contributions, a 
thorough numerical analysis must be implemented (it should be noted that this is the approach in Ref. 25] ). 

Since the disk and the cone have the same topology [ij] , the suppression of one degree of freedom in Eq. (ITTIl is a 
geometrical rather than a topological effect — it is due to a curvature singularity. 

Although the background addressed in this work is a cone, for which a is not greater than unity [cf. Eq. |^], it 



should be pointed out that in the context of topological defects in solids Mi?| there is room for a > 1 in the effective 
geometry corresponding to disclinations. In this case, as a increases from unity, Eii approaches Eqi and Oq, in Eq. 
(Pn|l decreases, suggesting that the classical limit in Eq. (|^ would be reached at lower temperatures (this point 
needs further investigation though). 
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